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Photoionization  of  Atoms  and  Ions: 

Application  of  Time-Dependent  Response 
Method  Within  the  Density  Functional  Theory 

I.  INTRODUCTION 

Accurate  calculation  of  photoionization  and  photoexcitation  cross- 
sections  of  atoms  and  ions  are  useful  in  a  variety  of  investigations  in 
plasma  physics  and  atomic  physics.  It  is  particularly  useful  in  the 
context  of  flashlamp  photopumping  schemes  for  x-ray  lasers.  Most  of  the 
existing  calculations  of  photoionization  cross-sections  were  done  using 
the  single  electron  or  the  independent  particle  model  (IPA).  In  this 
model,  the ‘energy-levels ,  and  vavefunctions  of  the  atom  or  ion  are  first 
calculated  using  the  Bartree-Fock  (HF)  method.  The  interaction  of  the 
incident  electromagnetic  radiation  with  the  atom  (or  ion)  is  treated  via 
the  first  order  perturbation  theory. 

Comparison  of  experimental  data  vith  the  IPA  calculations  shovs  that 
for  some  simple  systems  such  as  a  neutral  few  electron  atom  (Lithium,  for 
example),  there  is  qualitative  and  sometimes  quantitative  agreement. 
However,  for  many  electron  atoms  (Xenon,  for  example)  and  ions  vith  a 
large  number  of  bound  electrons,  substantial  discrepancies  are  found 
betveen  experimental  and  IPA-data.* 

In  our  work,  ve  used  the  time-dependent  linear  response  approximation 

vithin  the  framevork  of  the  relativistic  density  functional  method 
2  3  4 

(DFM)  ’  *  to  treat  the  problem  of  photoionization.  This  method 
incorporates  certain  advantages  over  the  HF-method.  The  HF-method  is  non¬ 
local  and  computationally  very  elaborate,  vhereas  in  the  density 
functional  method,  one  deals  vith  a  set  of  local  equations  only.  This 
leads  to  computational  simplicity.  On  the  other  hand,  it  is  veil  knovn 
from  extensive  application  of  the  density  functional  method,  that  fairly 
accurate  atomic  energy  levels,  vavefunctions  etc.  are  obtained.  The 
computational  simplicity  is  even  more  apparent  in  the  case  of  relativistic 
DFM  versus  relativistic  HF-methods.  In  the  DFM,  correlation  effects  of 
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the  bound  electrons  in  the  atoe  are  accounted  for  in  a  siaple  way  via  the 
correlation  potential.  The  Hartree-Fock  method,  on  the  other  hand,  does 
not  take  into  account  electron  correlation,  although  it  accounts  for  non¬ 
local  exchange  effects  appropriately. 

The  independent  particle  method  does  not  take  into  account  the 
polarization  effect  of  the  aton  brought  about  by  the  incident  time-varying 
radiation  field.  In  the  linear  response  method  within  the  density 
functional  method,  this  is  treated  adequately  -  as  will  be  seen  from 
comparison  with  the  experimental  data.  In  most  experimental  situations, 
the  incident  radiation  (from  synchrotron  sources  or  lasers)  have  field 
strengths  small  compared  to  the  atomic  field  strengths.  For  those 
experimental  conditions,  the  present  model  based  on  linear  response  is 
adequate  and  useful. 

Calculations  of  photoionization  and  photoexcitation  cross-sections  and 
rates  have  a  number  of  applications.  For  the  photopumping  scheme  for  x- 
ray  lasers,  these  processes  play  a  crucial  role  in  contributing  to  a 
population  inversion  of  excited  ionic  levels.  As  another  example, 
computation  of  opacities  of  plasmas  for  diagnostic  and  target  response 
effects  require  these  data  as  input.  In  order  to  model  the  radiation 
spectra  from  hot  plasmas  (via  detailed  configuration  rate  equation 
technique,  for  example),  the  photoionization  and  photoexcitation  data  are 
required  in  addition  to  other  bound-bound,  bound-free  and  free-free 
processes.  Acccurate  calculations  are  also  necessary  for  interpreting 
experimentally  available  data  on  cross-sections.  In  view  of  these 
different  applications,  there  is  a  need  for  realistic  modeling  of  these 
processes  in  order  to  generate  accurate  data  over  a  vide  range  of  photon 
energy  for  a  variety  of  atoms  and  ions.  The  present  model  provides  such  a 
tool  and  its  usefulness  will  be  discussed  in  subsequent  sections. 
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The  above  discussion  deals  with  the  calculation  of  photoionization  of 


isolated  atoms  or  ions  -  that  is,  without  the  effect  of  the  plasma 
environment.  These  effects  are  negligible  for  very  low  density  plasmas. 
In  many  experimental  situations,  however,  the  plasma  density  can  be  quite 
high  (for  example,  in  laser  produced  plasmas).  In  those  cases,  additional 
effects  due  to  screening  shifts  of  energy  levels  and  modification  of 
vavefunctions  and  potentials  for  the  ion  embedded  in  high  density  plasmas 
have  to  be  considered.  In  the  later  part  of  this  report,  we  will  present 
some  results  of  photoionization  cross-section  of  ions  in  a  dense  plasma 
medium  and  examine  the  modification  of  cross-sections  due  to  the  plasma 
environment. 

II.  THE  METHOD  OF  CALCULATION 

The  first  part  of  the  calculation  is  to  generate  the  energy-level 
spectrum  and  the  vavefunctions  of  the  particular  atom  or  ion  of  specific 
configuration.  This  is  done  by  using  the  local  density  functional  method. 
In  order  to  treat  many-electron  atoms  (vith  high  Z)  appropriately, 
relativistic  DFM  equations  vere  used.  In  this  method,  the  following  set 
of  equation  vere  solved  self-consistently: 


[c  a  .  p  ♦  c20  ♦  u  <r)]  ^  (r)  -  Ej  ^  (r)  (1) 

u  <r)  *  “  r  +  I  &d\'  +  [  p<?)  exc  <2> 

r  J  'r^r' 1  3p(r) L  xc  J 


and  p  (r)  -  I  f.|  ^(r)  j2 
i  1  1 


(3) 


In  the  above,  p(r)  is  the  electronic  charge  density  of  the  atom,  a's 
are  the  Dirac  matrices,  f^'s  are  the  integral  occupation  factors 
corresponding  to  the  number  of  electrons  in  each  state  ^(r)  vith 
corresponding  energy  eigenvalue  E^.  The  atomic  potential  u(r) 


ri-i 


contains,  in  addition  to  the  nuclear  and  the  electrostatic  Hartree 
term,  a  contribution  arising  from  the  electron  exchange  and 
correlation  effects.  Let  us  note  that  the  use  of  integer  occupation 
factors  f^'s  for  the  given  configuration  distinguishes  this  model  from 
the  "average  atom  model”  vhere  the  occupation  factors  are  taken  to  be 
those  given  by  the  statistical  Fermi  distribution  function. 

The  orbital  functions  are  four-component  spinors.  They  are 
split  into  major  and  minor  components: 


^(rn  r(A(r)/r)  i1  9j3jl(r) 


*2<r)J  l(B(r)/r)  i1  fljp^r) 


vhere  A  and  B  are  major  and  minor  components  of  the  radial  functions  and 
QJlm  and  are  tvo-component  Pauli  spinors  with  the  indicated  numbers. 
The  various  quantum  numbers  are  related  by 

1'  -  1  +  S,  J*  -  1  ♦  1/2,  S  •  1'  -  1/2  S,  K  -  -  S  (j  ♦  1/2);  S  »  ±  1  (5) 
The  differential  equations  for  A  and  B  (in  matrix  fora)  are: 


U) 


-  K/r 


(u  -  E  -  c  )/cs 


-  (u  -  B  c  )/cs 


In  Eq.  (2),  cxc  is  the  exchange-correlation  energy  of  the  electrons. 

3 

In  actual  calculation,  Gunnarsson-Lundquist  (G-L)  form  for  exchange- 
correlation  energy  and  potential  vas  used.  It  is  veil  knovn  that 
reliable  atonic  data  is  obtained  from  the  use  of  G-L  exchange- 
correlation.  Equations  (1)  -  (6)  are  solved  numerically  to  self- 

consistency  to  obtain  the  vavefunctions  ♦j's,  the  binding  energies  of 
each  orbital  E^,  the  atomic  charge  density  p(r)  and  the  self- 
consistent  potential  u(r). 
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Nov  consider  the  effect  of  an  incident  time-varying  radiation 
i  wt 

field  E  (t)  «  EQe  on  the  atom.  It  induces  a  time-dependent  atomic 
density  deviation,  6p(r, t),  causing  a  time-dependent  polarization 
effect.  For  the  linear  response  method  used  here,  it  is  convenient  to 
vork  vith  the  Fourier  transform: 

Mr,  t)  -  \-  J  Mr,«)e“iwt  d«  (7) 

The  net  induced  density  due  to  the  external  plus  the  induced  potential 
is 

Spind  "  J  Xfr.r' ,«)  [vext(r,«)  ♦  Vin<J  (r,«)j  dr',  (8) 

vhere  the  induced  potential  is  given  by 

sp<?:  •>  >vxc<p  <?»  (9) 

Vlmi<c'»>  *  J  j-;f  dr  *  - 

The  response  function  is  given  by 

X  (r,r',»)  -  E  fp+^r)  ♦i(r  )  G  (r,r  ,  E,  ♦  «) 

i  1  1  1  1  (10) 

♦  2  *i  ^(r)  (r  )  G*  (r,r  ,  Ei  -  ») 
i 

and  thus  involves  the  vavef unctions  and  energy  levels  of  the  atom. 
The  Green's  functions  are  solutions  of  the  inhomogeneous  Dirac 
equation 

(cm  .  p  +  c20  +  u  (r)  -  E)  G  (r,r  ,  E)  -  -  S(r  -  r  )  (11) 

In  actual  calculation,  angular  decomposition  of  the  Green's  function 
in  terms  of  spherical  harmonics  is  done  (Appendix  I)  and  the  radial 
part  is  treated  separately. 


The  frequency  dependent  polarizability  a(o>)  is  the  ratio  of  the 


induced  dipole  moment  to  the  xternal  field: 

<*(«)  -  -  J  Z  Sp  <r,«)  dr  (12) 

Note  that  ot(&>)  like  5p(r,w)  is  complex.  The  induced  density  deviation 
(ahd  also  the  corresponding  induced  potential)  can  have  a  phase 
difference  with  respect  to  that  of  the  applied  external  field.  Once 
a(»)  is  determined,  the  photoabsorption  cross-section  a( w)  of  the  atom 
is  obtained  from: 

e(«)  ,  15“  In  a  («).  (13) 

c 


III.  PARTIAL  CROSS-SECTION 

In  order  to  see  the  connection  with  the  IPA-model,  consider  the 
partial  cross-section  due  to  photoionization  from  a  specific  bound 
state  ^(r)  to  a  final  continuum  state  (r). 

The  initial  atomic  state  is  represented  as 


*i(r) 


(14) 


and  the  final  continuum  state  with  vavevector  it  and  energy  e  as 


*f(r)  -  Ait  Api1'  ?5il<£>  Y*,  (K)  Y^,  (r)  (15) 

t 

The  complex  coefficients  A^  s  are  found  by  requiring  <^(r)  to  behave 
asymptotically  as  an  incident  plane  vave  plus  a  spherical  wave.  Then 
the  partial  cross-section  is  shovn  to  be 


x  I  | A,  2| 
L'  A 


ani(«)  -  2  (21  ♦  1)  o-K  »  rT  a2 
<1  100!  1'0>2|  J  P£l,(r)  VSCF(r,u>)  URl(r)  dr | 2 


(16) 
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where  <1  100 |  1'0>  is  a  Clebsch-Gordon  coefficient. 

SCF 

In  (16),  V  (r,»)  is  a  frequency  dependent  complex  self- 

SCF 

consistent  potential.  Note  that,  if  V  (r,«)  is  replaced  by  the 
usual  dipole  moment  operator,  one  obtains  the  conventional  or 
independent  particle  approximation  (IPA)  result.  In  actual 

calculations,  both  bound  and  continuum  vavefunctions  are  generated 
numerically  using  the  Numerov  method  for  integrating  the  Dirac 
equation.  Let  us  also  note  that  the  real  and  imaginary  parts  of  the 
self-consistent  field  contributes  to  the  partial  cross-section  without 
interference.  Computations  were  performed  for  both  the  conventional 
independent  particle  model  and  the  time-dependent  linear  response  to 
density  functional  method  for  comparison  purposes. 

IV.  RESULTS 

In  Fig.  1,  results  of  computed  photoionization  cross-section  for 
neutral  xenon  is  plotted  as  a  function  of  photon  energy  near  the  Ad- 
threshold.  The  drastic  difference  between  the  results  obtained  from 
the  time-dependent  response  method  (curve  A)  and  from  conventional 
IPA-model  (curve  B)  is  clear.  The  IPA-model  does  not  reproduce  the 
experimental  values^  at  all  whereas  the  present  time-dependent  model 
agrees  very  well  with  the  -experimental  data  over  this  entire  range  of 
photon  energy  from  5  to  10  Ryd.,  including  the  peak  at  about  7  Ryd. 
In  that  energy  range,  the  IPA  cross-section  shows  a  rapid  decrease  - 
in  contrast  to  the  experimental  data. 

The  primary  physical  reason  for  this  difference  in  the  two  models 
arise  from  the  polarization  effect  of  the  atom  subjected  to  the 
incident  radiation.  In  the  case  of  xenon,  all  of  the  54  bound 


electrons  forming  the  atom  participate  in  the  polarization  process 
(the  outer  ones  contributing  most).  The  external  field  is  screened  in 
the  energy  range  5-6  Ryd.  (and  again  in  the  range  8.5  -  9.5  Ryd.),  and 
is  antiscreened  in  the  intermediate  6-8  Ryd.  range.  The  antiscreening 
effect  produces  a  stronger  effective  field  for  the  4d-electron  to 
photoionize,  thereby  enhancing  the  cross-sections  somewhat  in  the 
intermediate  energy  range  as  seen  in  Fig.  1. 

The  3s-partial  cross-section  for  Argon  is  shown  in  Fig.  2.  The 
experimental  data^  in  the  range  of  30  -  65  eV  are  depicted  by  circles 
with  the  error  bars.  The  conventional  IPA-calculation  (curve  A)  again 


does  not  shov  the  experimentally  seen  variation  at  all.  The  present 
time  dependent  response  method  reproduces  the  observed  variation 
including  the  experimental  Cooper  minimum  at  about  43  eV  and  is  in 


good  agreement  with  the  measured  cross-sections.  The  occurrence  of 

.1 


the  Cooper  minimum  is  known  to  be  due  to  the  vanishing  of  the  matrix 
element  between  bound  and  continuum  states  at  that  photon  energy  *  In 
the  present  time-dependent  polarization  model,  the  physical  reason  for 
the  minimum  is  that  the  induced  potential  almost  exactly  cancels  out 
the  external  field,  reducing  the  effective  field  to  almost  zero. 

The  total  cross-section  for  Argon  atom  as  a  function  of  the 


photon  energy  is- plotted  in  Fig.  3,  and  compared  with  experimental 
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data  .  Comparison  is  also  made  of  results  obtained  from  Hartree-Fock 

6 


(HF)  length  and  velocity  approximations  by  Kennedy  and  Hanson  .  We 
note  that  the  time-dependent  response  model  again  best  reproduces  the 
experimental  data.  The  HF-velocity  approximation  gives,  for  example, 
a  cross-section  twice  as  large  at  40  eV  whereas  the  HF-length  result 
is  six  times  larger  than  the  experimental  data  at  the  same  photon 
energy. 
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Good  agreement  with  available  experimental  data  for  Ar,  Ne,  Kr 
and  Xe  suggests  that  reliable  cross-sections  can  also  be  generated  for 
higher-Z  rare  gas  atoms  such  as  Radon  utilizing  the  time-dependent 
response  model.  For  Radon  (Z*86),  relativistic  effects  are 
significant.  The  Dirac  equation  approach  in  our  model  is  therefore 
suitable.  Ve  performed  self-consistent  calculation  for  Radon  and  the 
computed  total  cross-section  is  plotted  in  Fig.  4  for  the  photon 
energy  range  7-28  eV.  The  large  difference  betveen  the  results 
obtained  from  IPA-model  and  the  present  time-dependent  response  model 
is  clearly  seen  from  the  graph.  The  large  peaks  due  to  the  ionization 
from  ^3/2  *n  t*ie  IPA-model  are  masked  significantly  due  to 

polarization  effect  incorporated  in  the  time-dependent  response  model. 
No  experimental  data  is  available  for  the  case  of  Radon.  However,  in 
viev  of  the  good  agreement  betveen  experimental  data  and  the  results 
of  time-dependent  response  model,  it  is  expected  that  future 
experimental  measurements  for  Radon  in  this  photon  energy  range  will 
be  in  close  agreement  with  these  calculated  cross-sections. 

The  6s  partial  cross-section  for  Cesium  as  a  function  of  photon 

energy  in  the  range  of  4-8  eV  is  plotted  in  Fig.  5.  In  this 

particular  case,  there  are  significant  differences  in  the  two  sets  of 
9  10 

experimental  data  ’  of  the  two  groups,  presumably  due  to 
uncertainties  in  vapor  pressure  measurement.  Calculated  value  of 

partial  cross-sections  from  the  time-dpendent  model  is  in  fair 

agreement  with  one  set  of  observed  data  (Fig.  5).  The  effect  of  spin- 
orbit  interaction  is  known  to  be  significant  for  the  6s  partial  cross- 
section  of  Cesium.  A  more  detailed  theory  should  incorporate  this 
effect.  As  seen  from  Fig.  5,  the  IPA-model  shows  a  monotonic  decrease 
of  cross-section  in  this  energy  range  in  complete  disagreement  with 
the  experimental  data. 
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Ions  of  specific  configurations  such  as  Neon-like  Argon  are  of 
interest  for  scaling  to  higher  Z  neonlike  systems  for  x-ray  laser 
research.  No  experimental  data  is  available  for  these  ions  and  thus 
the  data  have  to  be  provided  by  theoretical  calculations.  The  fact 
that  the  time-dependent  response  method  is  successful  for  generating 
accurate  data  for  rare  gas  atoms  suggests  that  reliable  cross-sections 
can  be  generated  for  Ne-like  Ar  and  other  ions  of  interest.  Computed 
results  for  Ne-like  Argon  are  presented  in  Fig.  6  and  compared  with 
conventional  IPA-calculations.  The  differences  are  smaller  in  this 
case,  as  expected,  because  the  number  of  bound  electrons  for  these 
ions  is  only  ten.  For  these  highly  charged  ions,  these  electrons  are 
very  tightly  bound  and  thus  do  not  participate  as  effectively  in  the 
polarization  process. 

The  results  presented  above  are  for  single  atoms  or  ions  without 
the  effect  of  the  plasma  environment  -  as  appropriate  for  very  low 
density  plasmas.  For  high  density  plasmas,  however,  effects  due  to 
screening  shifts  of  energy  levels,  modification  of  wavefunctions  of 
bound  (particularly  the  upper  levels)  and  continuum  wavefunctions  as 
well  as  potentials  of  the  ion  embedded  in  the  plasma  have  to  be 
considered. 


For  proper  treatment  these  effects,  the  self-consistent  density 
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functional  method  (DFM)  at  finite  temperatures  should  be  used.  The 
application  of  this  method  requires  iterative  numerical  solution  of 
Schrodinger  equation  involving  the  complete  set  of  bound  and  continuum 
wavefunctions  for  the  raultielectron  ion  and  construction  of  effective 
potential  inclusive  of  plasma  screening  and  electronic  exchange 
correlation  effect  in  each  iteration.  For  our  present  purpose,  we 
adopted  the  following  simplified  approach  for  computational 
simplicity.  For  a  given  plasma  density  and  temperature,  we  represent 
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the  long  range  part  of  the  effective  potential  in  the  Debye-screened 
form.  The  inner  part  of  the  effective  potential  was  constructed  by 
numerical  itegration  of  the  Hartree  term  with  the  electron  density 
distribution  calculated  using  the  bound  state  wave functions. 

The  results  of  our  calculation  for  H-like  and  Li-like  carbon  ions 
are  shown  in  figures  7-11.  From  fig.  7 ,  we  see,  by  comparing  with 
isolated  atom  cross-sections,  that  the  threshold  for  photoionization 
from  ls-level  of  H-like  C  shifts  substantially  (from  18  eV  to  about 
15.2  eV)  for  Debye  length  Xp  -  2  a.u.  This  is  significant  in  the 
context  of  calculation  of  photoionization  rates,  which  is  calculated 
by  integrating  the  cross-section  weighted  by  the  electron  distribution 
function  (usually  the  Maxwell-Boltzmann  distribution)  over  the  entire 
energy  range.  The  calculated  rates  in  the  two  cases  would  therefore 
be  very  different  -  as  seen  from  fig.  7. 

In  fig.  8,  the  corresponding  results  of  partial  photoionization 
crosS-section  for  H-like  C  is  plotted.  Fig.  9  shows  the  3d-partial 
cross-sections  of  the  same  ion  at  two  different  plasma  conditions  (X^ 
-  5  and  2  a.u.)  and  compared  with  the  isolated  atom  cross-section 
(curve  A).  The  increasing  shift  of  the  ionization  threshold  with 
decreasing  Xp  (i.e.  with  increasing  density)  clear  from  the  figure  9. 

Ve  performed  similar  calculation  for  Li-like  carbon  as  well  (fig. 
10-11).  Since  there  are  three  bound  electrons  for  these  ions, 
iterative  calculations  mentioned  before  was  performed  for  the  bound 
states  and  the  continuum  vavefunction  was  computed  with  1-electron 
removed  (by  ionization)  from  the  given  state  2s  or  3s.  Comparison 


with  isolated  ion  results  shows  substantial  modification  of 


photoionization  cross-section  at  these  densities  and  temperatures  and 
significant  shifts  of  ionization  threshold  due  to  effect  of  the 
surrounding  plasma. 

V.  CONCLUSIONS 

It  is  demonstrated  that  the  time-dependent  linear  response  method 
within  the  framework  of  local  relativistic  density  functional  theory 
can  provide  reliable  atomic  data  for  various  atoms  and  ions  of 
experimental  interest.  This  model  is  particularly  useful  in  those 
situations  where  conventional  independent  particle  models  fail  to 
provide  accurate  data.  The  mechanism  of  time-dependent  polarization 
of  the  atom  is  seen  to  be  important  in  describing  the  observed 
results.  As  a  practical  point,  the  computer  code  based  on  the  time- 
dependent  model  is  fast  and  efficient,  capable  of  generating  a  large 
number  of  data  in  a  short  time  (for  example,  cross-sections  for  10 
photon  energies  for  a  medium-Z  atom  takes  about  3  minutes  of  c.p.u. 
time  on  a  Cray-XMP  computer).  The  present  method  is  capable  of 
treating  large  complex  atoms  with  high-Z  for  which  relativistic 
effects  are  important.  Let  us  point  out  that  if  the  applied  radiation 
field  strength  is  very  high,  so  that-  it  is  comparable  or  larger  than 
the  atomic  field  strength,  new  extensions  or  developments  of  the 
present  model  is  necessary  to  treat  those  conditions.  Full  numerical 
solution  of  time-dependent  density  functional  method  (beyond  the 
linear  response  approximation)  would  be  one  suitable  to  use  in  those 
cases.  Vork  in  this  new  direction  is  planned  for  the  future. 
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With  reference  to  the  calculations  for  the  different  plasma 
conditions,  let  us  point  out  that  the  Debye- screened  form  for  the  long 
range  part  of  the  potential  may  not  be  adequate  for  the  high-densi  ty 
plasmas.  For  more  realistic  calculation,  it  is  necessary  to  use  fully 
self-consistent  finite  temperature  density  functional  method  for  those 
plasma  conditions.  Also,  the  effect  of  exchange  and  correlation  for 
plasma  electrons  become  important  at  high  densities.  The  present 
approximate  scheme  is  used  at  present  for  computational  simplicity. 
Hovever,  the  results  shown  here  clearly  indicate  that  with  increasing 
plasma  density,  the  photoionization  cross-section  and  rates  of  the 
ions  forming  the  plasma  can  be  substantially  modified.  Accurate 
modeling  of  atomic  properties  for  dense  plasmas,  therefore,  requires 
that  the  effect  of  the  surrounding  plasma  should  be  properly  included 
in  the  calculation. 
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The  Green's  function  G  in  eq.  (11)  has  16  components,  vhich  are 


R 


K 


represented  in  matrix  form 


G(r,  r',E) 


G11(r, r' |  E)  G12(r,r'|  E)  ) 


l  G21(r, r' |  E)  G22(r,r' |  E)  ) 


The  angular  decompositon  of  various  terms  are 


<?,?•  |E)  .  aJUU>  g“  (r,r'  |E)  a‘u  (r') 


G12  <?,?■  | E)  .  jk  i1*1'  a,,.,<i)  G22(r,rt'f  |E)  at,,  (r-) 


‘jin"'  "jl" . .  *jl« 

G21  (?,?'  IE)  .  jf|t  I1'"1  Gn,„(0  G22(r,r’  |B)  a*la  ( i- ) 

G  22(rV  IE)  -  j{.  8n,.(0  G22,  (t.r«  |E)  a^,.  (t-) 


d-2) 


The  radial  part  G  (r,r'  E)  are  solutions  of  the  radial  inhomogeneous 
Dirac  equation 


Gjf(r,r' 


E,  -{ 


Ujl(r)  5jl  (C,)  '  r  <  r' 


Ujl(r)  ^1  '  (rr'Wjl)  C  >  r' 


(1-3) 


Vj ^  is  the  Vronskian 


[2  -  1  -  2  11 

j i ( r )  Uj^(r)  -  u j -^ ( r )  v j i ( r ) j  *  const. 


-1  -2 

(r)  and  (r)  are  major  and  minor  component  radial  functions  that  are 


real  and  regular  at  r  ■  o. 


-1  -2 

v  jl(r)  and  v  jj(r)  major  and  minor  ci'",r>onent 
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2 

radial  functions  which  (for  E  >  c  )  are  conplex  and  obey  outgoing  wave 

-1  -2 

boundary  conditions  at  r  »  •.  The  phase  for  u  and  u  ^  is  chosen  so 

2  -1  -2 

that  Wjl  is  real.  When  E  <  c  ,  Vjl(r)  and  Vj^(r)  are  real  and  decay 
exponentially  at  large  radii. 

With  the  above  representations,  the  polarizability  X(r,r'  u)  is  given 
by 


Xj  (r,r'  «)  -  I  fi 


U  1  2  1 

T  1  I 

1  I1  i2 

o 

o 

o 

.  1/2  j2  jr 

<2j2*l)  (211+1)  (212+1)  x 
2 


?"7'7E~;«T 
J2  2  1 


e^+  ») 


1 1  r  2  Bi<c»  "?2i2<r  ci  *  «■>]  *  [* 

i(r')  u  }  ,(r'  |et#) 

1  J2J‘2  1 

S  -S 

±  i  2  1  B.(r')  u  2  ,(r'|  e  ±  «)  / 

1  J2a2  1 

(rr')2|,  ] 

(1-4) 

vhen  r  <  r't  r  and  r'  are  interchanged  on  the  right  side  of  eq.  (1-4)  when 
r  >  r'.  The  index  i  stands  for  the  quantum  numbers  (n,  1^,  jj,  Sj)  of  a 
bound  state  and  f^  for  the  occupation  factors.  The  summation  is  over  all 
indices  except  I  and  over  both  +  u  and  -  «.  For  the  case  of  -  a,  the 


complex  conjugates  of  all  outgoing  waves  in  eq.  (1-4)  are  to  be  used. 
Angular  momentum  coupling  coefficients  are  expressed  in  terms  of  Vigner  3j 
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PHYS.  2,  178(1966).  ♦:  HARTREE-FOCK  (VELOCITY  APPROX.)  ANO 
0:  HARTREE-FOCK  (LENGTH  APPROX.)  BY  KENNEDY  &  HANSON. 

PHYS.  REV.  AS,  227(1972) . 


Figure  3 
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6S  PART IflL  CROSS-SECTION 


CES  UM 


PHOTON  ENERGY  (EV) 


SOLID  CURVE:  TIME-DEPENOENT  DENSITY  FUNCTIONAL  CALCULATION 
DASHED  CURVE:  INOEPENOENT  PARTICLE  MODEL. 

X:  EXPERIMENTAL  OATA  FROM  MOHLER  ANO  BOECKNER,  J.  RES. 

NBS  3,  303  (1929).  0:  EXPERIMENTAL  DATA  FROM  MARR  AND 

CREEK,  PROC.  R.  SOC.  A  304,  233  (1968). 
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CARBON  VI 


‘3D  LEVEL  CROSS  SECTIONS 


PHOTON  ENERGY  (eV) 

A:  ZERO  DENSITY,  B:  DEBYE  LENGTH  =  5  (AU) 
C:  DEBYE  LENGTH  =  2  (AU). 

Figure  9 
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PHOTOIONIZATION  CROSS-SECTION  (  x  1.0E-I8  CM**2) 


PHOTON  ENERGY  (eV) 

A:  ZERO  DENSITY.  B:  DEBYE  LENGTH  =5  (AU) 

Figure  11 
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